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ABSTRACT 



Aims. We present the first theoretical study on center-to-limb variation (CLV) properties and relative radius interpretation for narrow 
and broad-band filters, on the basis of a set of dynamic model atmospheres of C-rich AGB stars. We computed visibility profiles 
and the equivalent uniform disc radii (UD-radii) in order to investigate the dependence of these quantities upon the wavelength and 
pulsation phase. 

Methods. After an accurate morphological analysis of the visibility and intensity profiles determined in narrow and broad-band filter, 
we fitted our visibility profiles with a UD function simulating the observational approach. UD-radii have been computed using three 
different fitting-methods to investigate the influence of the sampling of the visibility profile: single point, two points and least square 
method. 

Results. The intensity and visibility profiles of models characterized by mass loss show a behaviour very different from a UD. We 
found that UD-radii are wavelength dependent and this dependence is stronger if mass loss is present. Strong opacity contributions 
from C 2 H 2 affect all radius measurements at 3 fim and in the N-band, resulting in higher values for the UD-radii. The predicted 
behaviour of UD-radii versus phase is complicated in the case of models with mass loss, while the radial changes are almost sinusoidal 
for the models without mass loss. Compared to the M-type stars, for the C-stars no windows for measuring the pure continuum are 
available. 

Key words. Stars: AGB and post-AGB - Stars: late-type - Stars: carbon - Stars: atmospheres - Techniques: interferometric 



1. Introduction 

The Asymptotic Giant Branch (AGB) is a late evolutionary stage 
of stars with masses less than about 8 M . The objects on the 
AGB are characterised by a degenerate C-O core and He/H- 
burning shells, a convective envelope and a very extended at- 
mosphere containing molecules and in many cases even dust 
grains. The atmospheres are affected by the pulsation of the inte- 
rior creating shocks in the outer layers. Due to the third dredge- 
up the AGB stars may have C/0> 1 dlben & Renzinil fl983) 
and their spectra are dominat ed by features of carbon species 
like C 2 , C 2 H 2 , C 3 , CN, HCN flGjagbgle t ynfT98llJovcdll998t 
Lancon et al.l l2000t iLoidl et al] 1200 U lYamamura & de Jongl 



needed for describing such objects i.e. dynamic model atmo- 
spheres. The status of modelling atmospher es of co o l AGB 
stars is summarised in different reviews by IWillsonl (|2000); 
IWoitkel d2003l) : iHofnetj d2005l 120071 l2009t) . Comparison of dy- 
na mic models with s pectroscopy for C-r i ch stars can be found 
in iHron et al.l dl998l): ILoidl et alj d200ll): iGautschv-Loidl et al. 
( 12004 henceforth GL04); INowotnv et al.l d2005allbl) : INowotnv 
d2005l) : lNowotnv et al.l d2009l) . 



2000). These are the classical Carbon stars. In such objects dust 



is mainly present as amorphous carbon grains. 

Studying stellar atmospheres is essential for a better com- 
prehension of this late stage of stellar evolution, to understand 
the complicated interaction of pulsation and atmospheric struc- 
ture, as well as the processes of dust formation and mass loss. 
The atmospheres of C-rich AGB stars with no p ronounced pul- 
sation can be described by hydrostatic models (J0rg ensen et alj 
2000; ILoidl et alj|200ll) . As the stars evolve the effective tem- 
peratures decrease and the atmospheres become more extended. 
At the same time the effects of time-dependent phenomena (dy- 
namic processes) become important and static models are not 
a good approximation anymore. More sophisticated tools are 



Because of their large atmospheric extension and their high 
brightness in the red and infrared, AGB stars are perfect can- 
didates for interferometric investigations. Moreover, since most 
observed AGB stars are very far away (at distance larger than 
100 pc), only the resolution reached by optical interferometry 
allows a study of the close circumstellar environment of these 
objects. 

The complex nature of their atmospheres opens the question 
of how to define a unique radius for these stars. In the literature 
several rad ius definitions were u sed, t hey are s umma rized in the 
reviews of lBaschek et al.l d!991l) and IScholzl d2003l henceforth 
S03). This is not a big problem when we deal with objects that 
have a less extended and almost hydrostatic atmosphere, because 
in most cases the different definitions of radius converge to a 
similar final result. However, one can find large differences when 
the envelope of the st ar is affected by the p henomena of pulsa- 
tion and mass loss. In Basche ket alj d!991l) the authors suggest 
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the use of the intensity radius that is "a monochromatic (or filter- 
integrated) quantity depending on wavelength or on the position 
and width of a specific filter". This radius is determined from 
the center-to-limb variation (CLV) as the point of inflection of 
the curve, in other words as the point where the second deriva- 
tive of the intensity profile is equal to zero. But to determine this 
value, a very accurate knowledge of the CLV has to be available. 
For AGB stars this is not always easy, since in the case of an 
evolved object the intensity profile can present a quite compli- 
cated behaviour and only few points are measured. A standard 
procedure in the case of a poorly measured CLV is to use analyt- 
ical functions like the uniform disc (UD), a fully darkened disc 
(FDD) or a gaussian to fit to the profile. The result of this fit is 
the so-called fitting radius (S03). 

Currently, some theoretical and observational studies about 
the different properties of the CLV and radius interpretations for 
M-type stars exist in the literature. They are summarised in the 
review of S03. According to this paper, the CLV shapes can be 
classified in: 

- small to moderate limb-darkening when the behaviour of the 
profile is well fitted by a UD; 

- the gaussian-type CLV is typical for extended atmospheres 
where there is a big difference in temperature between the 
cool upper atmosphere layers and the deep layers where the 
continuum is formed; 

- the CLV with tail or protrusion-type extension consists of 
two components, with a central part describable by the CLV 
of a near-continuum layer, and a tail shape given by the CLV 
of an outer shell; 

- the uncommon CLVs are profiles that show strange be- 
haviours and are characteristic for complex extended and 
cool atmospheres. 

The fitting radius is usually converted to a monochromatic 
optical depth radius R \, defined as the distance between the cen- 
tre of the object and the layer with t a = 1 at a given wave- 
length. But this conversion has to be done very carefully: even 
in the case of near-continuum bandpasses, spectral features can 
contaminate the measu re ment of the fitting rad ius as shown by 
Uacob & Scholzl d2002l) or lAringer etaT] d2008l) . and this will be 
not anymore a good approximation to the radius of the contin- 
uum layer. 

A filter radius, following the definition given by 
IScholz&Takedald 19871) 



Rfut = 



J R A I c J A dA 



(1) 



requires a good knowledge of the filter transmission curve and 
the filter itself has to be chosen in a careful way to avoid impure- 
filter-like effects (see Sect. 5 of S03 for more details). A very 
common definition of radius is the Rosseland radius that corre- 
sponds to the distance between the stellar centre and the layer 
with Rosseland otical depth tr oss = 1. As in the case of the 
monochromatic optical depth radius this quantity is not an ob- 
servable, it is model dependent, molecular contaminations are 
the main effect that has to be taken in account. 

A suitable window at 1.04 pm for determi ning the contin- 
uum r ad ius of M-type s ta rs was discussed by[Sc holz & Takeda 
d 19871). (Hofmann et alJ d!998l) . Uacob & Scholzl d2002l) and 
iTei et al.l d2003allbl) . The authors defined a narrow-bandpass lo- 
cated in this part of the spectrum, mostly free from contami- 
nations of molecules and lines. They indicate this radius as the 



best choice for studying the geometric pulsation. Unfortunately 
such a window does not exist in the case of C-rich AGB stars 
(cf. Sect. 0, and some model-assumptions have to be made for 
defining a continuum radius. On the other hand, available dy- 
namic atmospheric models for Carbon-stars dHofner et alf 2003) 
are more advanced than for M-type stars, the main reason being 
a better understanding of the dust formation. 

Following the same approach as done for M-stars by Scholz 
and collaborators, we started an investigation on the CLV char- 
acteristics fora set of dynamic models of pulsating C-stars 
dHron et alj|2008l) . 

In this work we present an overview on the main proper- 
ties of the intensity and visibility profiles in narrow and broad- 
band filters for dynamic model atmospheres with and without 
mass loss (Sect. [272] i. We use then the UD function to determine 
the UD-radius of the models, and investigate the dependence 
of the latter upon wavelength and pulsation phase. We discuss 
also a new definition for a continuum radius for the carbon stars 
(Sect. El). 



2. Dynamic models and interferometry 

2. 1 . Overview of the dynamic models 

The dynamic mode ls used in th is work are described i n detail in 
iHofner et al.l (I2003I) . GL04 and lNowotnv et all d2005al) . 

The models are spherically symmetric and characterised 
by time-dependent dynamics and frequency-dependent radia- 
tive transfer. The radiative energy transport in the atmosphere 
and wind is computed using opacity sampling in typically 50- 
60 wavelength points, assuming LTE, for both gas and dust. 
For the gas component, an energy equation is solved together 
with the dynamics, defining the gas temperature stratification 
(and allowing, in principle, for deviations from radiative equi- 
librium, e.g., in strong shocks). For the dust, on the other hand, 
the corresponding grain temperature is derived by assuming ra- 
diative equilibrium. Dust formation (in this sample only amor- 
phous ca rbon dust is considered) is treated by the "moment 
method" dGail & Sedlmavdll988tlGauger et al.|[l990l) . while the 
stellar pulsation is simulated by a piston at the inner bound- 
ary. Assuming a constant flux there this leads to a sinusoidal 
bolometric lightcurve which can be used to assign bolometric 
phases to each 'snapshot' of the model with = corresponding 
to phases of maximum light. 

In Table Q] the main parameters of the models used here are 
summarised. Luminosity L*, effective temperature and C/O 
ratio refer to the initial hydrostatic model used to start the dy- 
namic computation. This model c an be compared with c lassic 
model atmospheres like MARCS dGustafsson et alJl2008T) . The 
stellar radius of the initial hydrostatic model is obtained from 
L* and and it corresponds to the Rosseland radius. The piston 
velocity amplitude Au p and the period P are also input param- 
eters and rule the pulsation, while the mass loss rate M and the 
mean degree of condensation (/ c ) are a result of the dynamic 
computations. The free parameter /l is introduced to adjust the 
amplitude of the luminosity variation at the inner boundary in a 
way that L m {t) °c f L ■ R^t) (GL04). 

The models used in this work were chosen with the purpose 
to investigate the effects of the different parameters on radius 
measurements: mainly mass loss but also effective temperature, 
C/O and piston amplitude. The sample can be divided into two 
main groups, keeping this purpose in mind: models with and 
without mass loss (in TableQ]series D and N, respectively). 
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Table 1. Parameters of the dynamic model atmospheres used in this work (see also GL04). Luminosity (L»), temperature (r*), 
Rosseland radius (/?*) and C/O refered to the initial hydrostatic model, amplitude of the piston velocity (Am p ), period (P) and /l are 
input parameters for the dynamic calculations, while mass loss rate (M) and the mean degree of condensation (/ c ) are a result. All 
the initial models have a mass M* = 1 M G and solar metallicity. Models in bold are the reference chosen for this work. See also text 
(Sect.©. 



Model 




T, 


R, 


C/O 


Au p 


P 


,/l 


AMboi 


</c> 


M 








[K] 


[Re] 




[kms- 1 ] 


[d] 




[mag] 




[Mo 


yr" 1 ] 


Dl 


10000 


2600 


492 


1.4 


4.0 


490 


2.0 


0.86 


0.28 


4.3- 


10 6 


D2 


10000 


2600 


492 


1.4 


4.0 


525 


1.0 


0.42 


0.37 


5.9- 


io- 6 


D3 


7000 


2600 


412 


1.4 


6.0 


490 


1.5 


1.07 


0.40 


2.5 • 


10~ 6 


D4 


7000 


2800 


355 


1.4 


4.0 


390 


1.0 


0.42 


0.28 


2.4- 


10~ 6 


D5 


7000 


2800 


355 


1.4 


5.0 


390 


1.0 


0.53 


0.33 


3.5- 


10~ 6 


Nl 


5200 


3200 


234 


1.1 


2.0 


295 


1.0 


0.13 








N2 


7000 


3000 


309 


1.1 


2.0 


390 


1.0 


0.23 








N3 


7000 


3000 


309 


1.4 


2.0 


390 


1.0 


0.24 








N4 


7000 


3000 


309 


1.4 


4.0 


390 


1.0 


0.48 








N5 


7000 


2800 


355 


1.1 


4.0 


390 


1.0 


0.42 









It was already mentioned that the pulsation of the inner lay- 
ers can influence the outer atmosphere of the star. The propaga- 
tion of the shock waves causes a levitation of the outer layers 
and in certain cases formation of dust grains can take place. The 
radiation pressure on the dust grains may drive a stellar wind 
characterized by terminal velocities and mass loss rates M. 
In Fig.[TJan example of the radial structure of the dynamic model 
atmosphere Dl (scaled to R+) is plotted for different phases. For 
comparison, the area covered by the model Nl is also shown 
grey-shaded. In particular we show the effect of the dust-driven 
wind that appears on the model structure Dl. The different at- 
mospheric extensions of the two models are easily noted. 

In the series D of Table[TJtwo models, Dl and D2, differ from 
the others by a higher luminosity (L*). The only distinction be- 
tween these two models are different values of the period P and 
the parameter /l. The effect on the structure of the atmosphere 
is that dust shells observed at few R/R+ are more pronounced 
in the case of D2. The models D4 and D5 have the same lumi- 
nosity, temperature and period but different piston velocities and 
different resulting degrees of dust condensation. In the structure 
of D4, compared with D5, there is a smoother transition from 
the region dominated by pulsa tion to the one domina ted by the 
dust-driven wind (see Fig. 2 of lNowotnv et all2005al) . 

In the N set of models we can distinguish Nl by its higher 
temperature and lower luminosity. All the other models have the 
same luminosity but differ by other parameters. N3 and N4 are 
an example for studying the effect of different piston velocities 
on models without mass loss. The atmospheric stucture of N4 is 
more extended but the outer layers are not yet cool and dense 
enough for dust formation. N2 and N3 have a different C/O ratio 
that causes a difference in the spectral features but also in the 
estimation of the diameter at certain wavelengths (Sect. [5]). 

We will use model Dl as reference among the ones with 
mass loss. This model was compared in GL04 with observed 
spectra of the star S Cep, it was shown that different phases of D 1 
can reproduce the S Cep spectra shortward of « 4 /mi. The same 
model is in qualitative agreement with observ ed CO and CN 
line pr ofiles variations of S Cep and other Miras (Now otnv et alj 
2005a b). Concerning models without mass loss we use model 
Nl as a reference, which gives the best fit of spectra and colour 
measurements for TX Psc in GL04. Each one of the reference 
models represents also the most extreme case in its sub-group. 



2.2. Synthetic intensity and visibility profiles 

To compute intensity and visibility profiles on the basis of the 
models we proceeded in the following way. 

For a given atmospheric structure (where temperature and 
density stratifications both are taken directly from the dynami- 
cal simulation) opacities were computed using the COMA code, 
originally developed by Aringer (2000): all important species, 
i.e. CO, CH, C 2 , C 3 , C 2 H 2 , HCN and CN, are included. The data 
so urces used for the continu ous opacity are the ones presented 
in lLederer & Aringer! (12009). Voigt profiles are used for atomic 
lines, while Doppler profiles are computed for the molecules. 
Conditions of LTE and a microturbulence value of £ = 
2.5 km s are assu med. All references to the molecular data are 
listed in Table 2 of lLederer & Aring er (2009 ). The sourc e for the 
amorphous carbon opacity is iRouleau & Martini d 199 lb. These 
input data are in agreement with the previous works dNowotnvl 
2005; iLederer & Aringer! 120091; lAringer et all 120091) . and with 
the opacities used for the computation of the dynamic models. 
Information on recent up dates on the COMA code and opac- 
ity da t a can be found in iNowotnvl (120051). ILederer & Ar inger 
(l2009h . lAringer et alJd2009l) and lNowotnv et al.ld2009h . Theres- 
olution of the opacity sampling adopted in our computations is 
18 000. 

The resulting opacities are used as input for a spherical radia- 
tive transfer code that computes the spectrum at a selected reso- 
lution and for a chosen wavelength range. In addition, the code 
computes spatial intensity profiles for every frequency point of 
the calculation. Velocity effects were not taken into account in 
the radiative transfer. 

The monochromatic spatial intensity profiles were subse- 
quently convolved with different filter curves (narrow and broad- 
band) in the near- to mid-IR range to obtain averaged inten- 
sity profiles. The narrow-band filters are rectangular (transmis- 
sion 1), while the broad-b and filters are the standard filters 
from Bessel & Brettj dl988l) . The visibility profiles are the two- 
dimensional Fourier transform of the intensity distribution of the 
object, but under the condition of spherical symmetry (as for our 
ID models), this reduces to the more simple Hankel transform 
of the intensity profiles (e.g. Bracewell 1965). 

The intensity and visibility profiles in the case of the broad- 
band filters were computed with the following procedure. Each 
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1 10 
R/R, 



Fig. 1. Radial structure of the dynamic model Dl. Gas tempera- 
ture (upper panel) and density (lower panel) are plotted for dif- 
ferent phases during the pulsation period. For comparison the 
grey-shaded area marks the range where the structure of the 
model Nl at different phases can be found. More details about 
the structure of the m odels can be found in the text and in 
iNowotny et all (l2005allbh . 



Table 2. Definitions of the narrow and broad band filters. Listed 
are the names of the bands, characteristic features (for narrow- 
bands), central wavelengths, bandwidths and the morphological 
classes (see Sect. 01. In bold are indicated the filters chosen to 
represent the respective morphological class. The designation 
"cont" means that the spectrum includes no strong molecular 
absorption, but this does not necessarily represent the real con- 
tinuum. 



Filter 




A c 


A/1 


Class 






[/JITl] 






1 01 

l.Ul 


C2, CN 


1 D1 




j 


1 OQ 




1 OQ 

1 .U7 


0.02 


j 


1.11 




1.11 


0.02 


j 


1 51 


CN+ 


1 51 

1 .J 1 


0.02 


j 


1.53 




1.53 


0.02 


I 


1.575 


CN 


1.575 


0.05 


I 


1.625 


CN, CO 


1.625 


0.05 


I 


1.675 


"cont" 


1.675 


0.05 


I 


1.775 


^2 


1.775 


0.05 


I 


2.175 


CN 


2.175 


0.05 


I 


2.375 


CO 


2.375 


0.05 


I 


3.175 


C 2 H 2 , HCN, CN 


3.175 


0.05 


III 


3.525 


"cont" 


3.525 


0.05 


11 


3.775 


C2H2 


3.775 


0.05 


II 


3.825 


C2H2 


3.825 


0.05 


II 


3.975 


"cont" 


3.975 


0.05 


II 


9.975 


"cont" 


9.975 


0.05 


II 


11.025 


"cont" 


11.025 


0.05 


II 


12.025 


C2H2 


12.025 


0.05 


111 


12.475 


C 2 H 2 


12.475 


0.05 


111 


12.775 


C2H2 


12.775 


0.05 


III 


J 




1.26 


0.44 




H 




1.65 


0.38 




K 




2.21 


0.54 




L' 




3.81 


0.74 





In the following two sections we will discuss the character- 
istics of the intensity and visibility profiles for different mod- 
els and the effects when using different filters. In the case of 
broad-band filters also the so called b and width smearing effect 
( Verhoelst 2005; K ervella et alj|2003l) will be discussed. 



broad-band was split into a set of 10 narrow-band subfilter£] and 
the squared broad-band visibility was computed as 

y 2 — _! cy\ 

'broad VCC2f^ ' ' 



where the sum is calculated over all the narrow subfilters. S \ is 
the transfer function of the broad-band filter within the corre- 
sponding narrow -band subfilter, F, is the flux integrated over the 
corresponding narrow-band subfilter, and V,- is the visibility of 
the narrow-band subfilter computed from the average intensity 
profile for the corresponding filters. 



1 In principle this should be done monochromatically but for compu- 
tational reasons the usual number of narrow-band subfilters was limited 
to 10. However, we compared the resulting broad-band filters computed 
with a set of 10 narrow-band subfilters with the one computed with a 
set of 100 and 200 subfilters without finding any significant difference 
in the resulting broad-band visibility profile. We want to stress also that 
the maximum resolution for the computation of the broad-band filter is 
given by the opacity-sampling resolution. 



3. Synthetic profiles for narrow-band filters 

Aiming to study effects of different molecular features charac- 
teristic for C-star spectra and somewha t following the approach 
sketched in iBessel etaTI dl989, 1996), 21 narrow-band filters 
have been defined in the near- to mid-IR range with a "resolu- 
tion" ranging from 50 to 200. Figure|2]shows the location of the 
filters on top of a typical synthetic (continuum normalised) spec- 
trum of a C-type star for illustration. Table [2] gives an overview 
of our complete set of filters with central wavelengths, band 
widths, and related spectral features. A few of these filters are 
labeled "cont" meaning that the contamination with molecular 
features is comparably low in this spectral region. It may be pos- 
sible at extremely high spectral resolutions to find very narrow 
windows in the spectrum of a C-star where no molecular lines 
are located (e.g. Fig. 6 in Nowotny et al. 2005b). However, in 
general and with resolutions available nowadays the observed 
spectra of evolved C-type stars in the visual and IR are crowded 
with absorption features. Due to such a blending the continuum 
level F/F Q -\ is never reached as in Fig. [2] This effect is more 
and more pronounced with lowe r effective temperatur es and de- 
creasing spectral resolutions (cf. Now otny et al.l l2009). 
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J I I L 



.8 1 



10 





1 n 1 




111 




I . I I 

1 S 1 
i .J i 




1 ^ 7 




1 




1 R9S 




1 R7S 




1 77S 




9 1 7S 

Z. • I / U 








1 7S 




J. JZ J 




3.775 




3.825 




3.975 




9.975 




1 1.025 




12.025 




12.475 




12.775 


18 



[yum] 



Fig. 2. Synthetic spectrum (R = 360) calculated on the basis of the hydrostatic intial atmospheric structure of model D 1 . All nece ssary 
opacity sources (continuous, atoms, molecules, dust) are included for the spectral synthesis, (see Fig. 4 of lAringer et alj d2009l) for 
the individual contributions of the different molecular species). The flux is normalised to a calculation for which only the continuous 
opacities were taken into account. Overplotted are the narrow -band filters as defined in Table 2, for which the central wavelengths 
in [/urn] are given in the legend. In addition, the ranges of the broad-band filters are marked grey-shaded.. 



As a first step we study the morphology of different inten- 
sity and visibility profiles, then we proceed with an investigation 
of the main opacity contributors to the features that appear in 
the profiles. For this purpose we compute for the two reference 
models three different sets of intensity and visibility profiles: 
(i) profiles including all necessary opacity sources i.e. continu- 
ous, atomic, molecular and dust opacities; (ii) profiles where the 
opacities of dust were neglected, to analyse only molecular and 
atomic contributions, and (iii) profiles for which only the con- 
tinuous gas opacity is considered, to get theoretical continuum 
profiles. 

For the narrow-band filters, the intensity profiles of models 
with mass loss can be divided into three morphological classes 
listed in the last column of Table [2] It turns out that this classi- 
fication is related mainly to the wavelength region with an ex- 
ception for the filter 3.175 /mi which belongs to the third class 
(that contains filters centered around 12 /mi) instead of the sec- 
ond (that contains all the other filters in the 3 //m region plus the 
9 and 1 1 /mi filters). In Fig.[3]the intensity and visibility profiles 
for the representative filters of each class are illustrated. Models 
with mass loss (left panels) are more extended than models with- 
out mass loss (right), and their intensity profiles fit into the cate- 
gory of "uncommon profiles" of S03. 

The first class of profiles contains the filters in the near IR 
from 1.01 jt/m to 2.375 /mi. It can be well represented by filter 
1.53 (C2H2, HCN, CN). Its morphology is characterised first 
by a limb darkened disc, followed by a bump due to molecular 
contributions, and a peak indicating the inner radius of the dust 
shell. At 1.53 /mi a feature due to the molecular species of C2H2 
and HCN is predicted by some of the cool dynamic models with 
wind (GL04). This feature was observed in the spectra of cool 
C-rich Miras such as RLep dLancon et alj|2000l) . VCyg, cya41 



(faint Cygnus fiel d star) and s ome other faint high Galactic lati- 
tude carbon stars (Joyce 1998). 

The third class of profiles is represented by filter 3.175 
(C2H2, HCN), which is placed right in the center of an ex- 
tremely strong absorption feature characteristic for C-star spec- 
tra (Fig. [2j. Observations in this spectral region are a very pow- 
erful tool for understanding the upper photosphere. Here, the in- 
ner central part of the intensity profile can be described by an 
extended plateau due to strong molecular opacity. All the filters 
defined around 12 /mi - where C2H2 has an enhanced absorption 
- belong also to this class. 

The second class of profiles is intermediate between the two 
behaviours illustrated before with an indication of an outer dust 
shell around 8 AU. It is represented by filter 9.975 ("cont"), and 
it covers the behaviour of the filter at 1 1 .025 /mi plus those de- 
fined in the 3 /mi region (except 3.175 /mi). 

In the Fourier space the three representative visibility pro- 
files are also markedly different (lower left panel of Fig. [3). The 
UD function, computed by fitting the point at visibility^ V = 0.3 
of the 3.175 /mi profile (for details see Sect. [Sj^ is shown for 
comparison in Fig. [3] The visibility profiles from the model and 
the UD are different at visibility < 0.2, in this region the model 
profile is less steep and has a less pronounced second lobe com- 
pared with the UD. 

In the right panel of Fig.[3]intensity and visibility profiles for 
the phase of maximum luminosity of the reference model with- 
out mass loss are represented. As expected from the model struc- 
tures, the intensity profiles are more compact than the ones for 
the models with mass loss. In the intensity profile corresponding 
to the filter 9.975, around 1.3 AU, a peak due to the molecu- 
lar opacity of C3 is visible. J0rge nsen et alj d2000l) showed that 

2 All the visibility values V have to be interpreted as normalized vis- 
ibilities V/V . 
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Fig. 3. Intensity and visibility profiles for the narrow-band filters at the maximum phase of the mass-losing model Dl on the left 
side, and the model without mass loss Nl on the right side. All the opacities are included in the computations. Profiles for three 
representative filters are shown together with the equivalent UD function for the 3.175 yum filter. The UD is fitted to the visibility 
point V = 0.3 marked with the diamond (details in Sect.|5}. 



C3 is particularly sensitive to the temperature-pressure structure 
and, hence, tends to be found in narrower regions than other 
molecules. The peak visible at 3 AU of the profile of the model 
Dl is due mainly to C3 but also to other molecular and dust opac- 
ity contributions. These 9.975 fim peaks in the intensity profiles 
mean a marginally "ring"-type appearence of the disc. The shape 
of the intensity profiles for the models without mass loss can be 
classified as having small to moderate limb-darkening (S03). In 
the Fourier space, the visibility profiles of the three filters for the 
series N look very similar except for tiny differences especially 
in the second lobe. In all the cases they can be well approximated 
by a UD in the first lobe. 

Plots of theoretical continuum profiles for both reference 
models are very similar to UDs. The classification of filters is 
valid also for all the other models of our sample listed in TableQ] 



4. Synthetic profiles for broad-band filters 

We investigated also intensity and visibility profiles for the 
broad-band filters J, H, K, L', because several observations of 
C-stars exist in these filters. 

In Fig. |4] the resulting broad-band visibility profiles com- 
puted for the maximum phase of the reference models are shown 
together with the corresponding UD visibility, which was cal- 



culated with the one-point-fit method at V = 0.3. In the case 
of model Dl (with mass loss) the broad-bands profiles show a 
behaviour completely different from the UD. The slopes of the 
model profiles are always less steep than the corresponding UD 
at V < 0.2. The synthetic visibility profile shows a characteristic 
tail-shape instead of the well pronounced second lobe that com- 
pare in the UD. The broad-band visibilities for the model Nl are 
well represented by a UD and, like in the case of the narrow- 
bands, the major (but small) differences occur in the second lobe 
of the profiles. 

When a comparison to observations is performed in a broad- 
band filter particular attention has to be paid to the smearing 
effect (also kn own as bandwidth smearing) which affects the 
measurements dDavis et alj2000llKervella et alj2003tlVerhoelstl 
2005). This is nothing more than the chromatic aberration due 
to the fact that observations are never monochromatic but in- 
tegrated over a certain wavelength range defined by the broad- 
band filter curve. Furthermore, each physical baseline B, corre- 
sponds to a different spatial frequency u with u = j . A broad- 
band visibility, calculated at a given projected baseline, will re- 
sult in a superposition of all the monochromatic profiles over the 
width of the band. 

According to iKervella et ail d2003l) and IWittkowski et al.1 

(2004), the bandwidth smearing is stronger at low visibilities 
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Fig. 4. Visibility profiles computed in broad-band filters for the maximum phase of the model Dl with mass loss on the left side, 
and the model Nl without mass loss on the right side. The solid line corresponds to the model while the dashed line is the equivalent 
UD computed fitting one point at visibility V = 0.3. 



close to the first zero of the visibility profile. We will demon- 
strate that also the shape of the profile plays an important role in 
the smearing effect. 

We investigated this effect with the help of our reference 
models in the J, H, K and L' filters. As explained in Sect. 12.21 
the broad-band profile was computed splitting the broad-band in 
10 narrow-band sub filters and then averaging the resulting pro- 
files. To take in account the bandwidth smearing, we split again 
the broad-band in the same set of 10 narrow-band subfilters. The 
spatial frequencies are then converted from AU~ l to baselines 
in meters, fixing the distance for an hypothetical object (500 
pc) and using as corresponding "monochromatic wavelength" 
the respective central wavelength of the narrow -band filter. Then 
the resulting broad-band profile was obtained using the above 
Eq. ©. 

Figure [5] illustrates the bandwidth smearing effect for the 
Kband for the two reference model atmospheres (Dl on the 
left and Nl on the right). It is clearly visible in the upper pan- 
els that the profiles of the narrow-band subfilters reach the first 
zero at different baselines. The bandwidth smearing affects the 
br oad-band profile in a way that it never reaches zero, as noted 
in iKervella et al.l d2003). We compared the resulting visibility 
curve with a corresponding one obtained without taking into ac- 
count the bandwidth effect. The result of this experiment can be 
seen in the lower panels of Fig.|5]where we plot the difference (in 



percentage) between the visibility with and without taking into 
account the bandwidth smearing (BW and NBW respectively). 
For both types of models - with and without mass loss - the 
smearing gives a relative difference in the profiles that reaches 
at certain baselines more than 70%. For model Nl the major 
difference occurs near the first and the s econd zero and it is 
negligible elsewhere (a s already shown in IKervella et al.| [2003; 
Wittkowski et al. 2004), while in the case of the model Dl the 
behaviour of the relative difference is less obvious. This result is 
confirmed also for the other broad-band profiles. However, it ob- 
viously depends on the shape of the visibility profile, and hence 
on the model. 

We conclude that in all the cases it is important to compute 
accurate profiles taking into account the bandwidth smearing. 



5. Uniform disc (UD) diameters 

For all the models of the sample and all the considered filters we 
computed the equivalent UD-radii using three different methods: 
(i) a single point fit at visibility V = 0.3, (ii) a two-points fit at 
visibilities V = 0.1 and 0.4, and (iii) finally a non linear least 
square fit using a Levenberg-Marquardt algorithm method, and 
taking into account all the po ints with vi s ibility V > 0.1. As al- 
ready stressed in Sect. 5 of lire land et al.l d2004al) . none of these 
procedures are fully consistent because the UD does not repre- 




Fig. 5. Representation of the bandwidth smearing effect for the maximum of phase of the model Dl (left panels) and Nl (right 
panels) in the K band. In the upper panels visibility profiles as function of baseline for the individual narrow -band subfilters (dotted 
lines) and for the resulting broad-band filter (solid line) - computed taking into account the bandwidth smearing - are plotted. In 
the inset, a zoom of the region around the first zero is represented. The lower panels show the visibility profiles computed with 
bandwidth smearing (BW, solid line), and profiles computed neglecting the effect (NBW, dotted line). The dashed lines present the 
difference between the two profiles (scale on the right y-axis). 



sent the real shape of the profile (especially in the case of models 
with mass loss). However the UD is often chosen to interpret in- 
terferometric data with a small sampling of the visibility curve. 
Thus, the three different methods are used in this work to inves- 
tigate the influence of the sampling of the visibility profile, in 
particular if only a few visibility measurements are available. 
We use the visibility V = 0.3 for the one-point-fit following 
llreland et alj d2004allbl) who identified this point as a good esti- 
mator for the size o f the object. While thi s is true for the M-type 
models presented in llreland et al.l d2004allbl) . it is not obvious for 
our C-rich models with mass loss, which are very far from a UD 
behaviour (Sect. [3). However keeping in mind that a UD is just a 
first approximation we used this value also for our computations 
of the one-point fit. 

The aim of this section is to study the behaviour of the UD- 
radius versus wavelength and versus phase, as well as the depen- 
dence on the different parameters of the models. 

As we pointed out previously, while for the M-type stars 
the continuum radius can be measured around 1.04 fim (e.g. 
Uacob & Scholzl |2002|) . our calculations show that the C-stars 
spectra do not provide any suitable window for this determina- 
tion. However, the continuum radius is a relevant quantity for 



modeling of the stellar interior and pulsation. Therefore we com- 
puted for our sample of dynamic model atmospheres a set of 
profiles taking into account only the continuous gas opacity. The 
resulting profiles correspond to the pure theoretical continuum 
absorption of the models. An equivalent UD-continuum radius 
has been computed for this theoretical "pure continuum" pro- 
files. This allows us to study how the UD-continuum radius de- 
pends upon wavelength and phase and may give a possibility to 
convert observed radii to the continuum radius. 

5.1. Uniform disc radii as a function of wavelength 

The UD-radii resulting from least square fits at the maximum 
phase of all the models are plotted versus wavelength in Fig. [6] 
They are scaled to to minimize effects of the initial stellar 
parameters and enhance effects due to dynamics and dust for- 
mation. The UD-continuum radii are represented with the cross 
symbols for Dl (upper panel) and Nl (lower panel). The ma- 
jor difference is found between models with and without mass 
loss. In both cases the UD-radius generally increases with wave- 
length, but this behaviour is more pronounced for mass-losing 
models. At 3.175 /mi there is a notable "jump" (value signifi- 
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cantly higher than the surrounding) of the UD-radius due to the 
high C2H2 opacity. The UD-continuum radius is mostly indepen- 
dent of wavelength and thus it can be adopted as a good refer- 
ence radius (although not being observable for C-stars!). In fact, 
a small increase of its size can hardly be distinguished at longer 
wavelengths (see Fig. |6j, and at 1.6 /mi is visible the minimum 
opacity of H . The ratio between the apparent size of the contin- 
uum radius at 1 1 and and 2.2 /jm is 1 . 14 in the case of the model 
with mass loss Dl and 1.06 for the model without mass loss Nl. 
This change in apparent size can b e explained with the elec tron- 
hydrogen continuum, according to lTatebe & Townesl d2006l) . All 
the models show a value of UD-continuum radius lower than the 
stellar radius This is a consequence of the definition of R+. 
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Fig. 6. UD-radii versus wavelength for the maximum phase of 
different models with (upper panel) and without (lower panel) 
mass loss. The cross symbol represents the UD-continuum ra- 
dius computed for the two reference models (Dl, Nl). The 
method used for the UD determination is a least-square fit us- 
ing all the points of the synthetic profile with visibility > 0.1. 

Comparing the effect of the different model parameters we 
can say that model extensions increase with decreasing effec- 
tive temperature and higher mass loss values. The models Dl 
and D2 for example are distinguished by a higher mass loss rate 
for D2 because of the different period and /l value. The result- 
ing UD-radii for D2 are systematically higher than the one for 
D 1 . D5 has a higher piston velocity, and consequently a higher 
mass loss, compared to D4. The UD-radii obtained are larger. A 
comparison of D1-D2 or D4-D5 in Fig. [6] shows that the effect 
of the mass loss rate on the UD-radius is larger beyond 3 /mi. 



The model D3 has UD-radii systematically larger than the one 
of the model D4 in the region between 1 and 2 /zm and at 9 /mi. 
At 3.175 //m the trend is inverted while for the other filters at 3 
/mi and for the longer wavelength the resulting UD-radii is the 
same for the two models. D3 and D4 differ in the temperature 
(2 600 K for D3 and 2 800 K for D4), piston velocity, period and 
/l while the resulting mass loss rate is quite similar. 

On the other hand, in absence of dusty winds the models 
with the same parameters except for C/O (N2 versus N3) have 
very similar radii. The difference is barely detectable at 3 /mi and 
at longer wavelengths. If, besides different C/O ratios, also the 
temperature of the model changes (e.g. as between N4 and N5) 
the resulting UD-radii show more pronounced differences. The 
UD-radii of Nl overlap always with the one of N2, even if the 
models have different temperature (higher for Nl), luminosity 
(higher for N2) and period (longer for N2). N3 and N4 have all 
the same parameters except the piston velocity and the resulting 
UD-radii is larger for the second model, which has higher value 
of Aiip. Model N4 is also the most extended among its subgroup. 

5.2. Separating models observationally 

The dynamic model atmospheres provide, compared with clas- 
sical hydrostatic models, a self-consistent treatment of time- 
dependent phenomena (shock propagation, stellar wind, dust 
formation). It is thus of course highly interesting if observed UD- 
radii at different wavelengths can be used to distinguish between 
the different models and then to assign a given model to a given 
star. 

If the model UD-radii are not scaled to their respective 
the D models differ by 30% to 50%. These differences should be 
detectable given typical distance uncertainties of 20%. However, 
as dicussed above, this would mostly separate effects due to the 
fundamental stellar parameters (initial radius, luminosity) and 
allow to discriminate between N and D models but dynamic ef- 
fects within each series are only of the order of 10% to 25%. This 
can be seen from Fig. [6] or if one scales each model to e.g. its 
UD-radius in the K-band. The latter approach would be the most 
sensitive one as it eliminates errors in the distance (provided the 
measurements were done at the same pulsational phase, better at 
the same time!). Assuming that individual UD-radii can be de- 
termined with an accuracy of 5% to 10%, the dynamic effects 
could be (barely) distinguishable for the D models. The most 
sensitive wavelengths in this respect would be filters in the K- 
and L-band. The A^-band shows a similar sensitivity to dynamic 
effects though at a reduced observational accuracy. For the N 
models these effects are less then 10%, with the largest differ- 
ence again beyond the /T-band. 

Nevertheless, the best strategy for choosing the right dy- 
namic model is to combine spectroscopic and interferometric 
observations. 

Finally it should be noted, that for a larger grid of models 
(e.g. iMattsson et alj|2009h a wider separation of models might 
show up. This will however be the subject of a future paper. 

5.3. UD-radii as a function of time 

The behaviour of the UD-radius versus phase predicted by the 
models is rather complex, especially in the case of mass-losing 
models. In Fig.Qthe UD-radius at a given wavelength/filter ver- 
sus phase is plotted for the two reference models. The presence 
of molecular and dust opacities increases the UD-radius strongly 
compared to the UD-continuum radius. In both types of models 
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Fig. 7. UD-radii versus phase are shown for three different narrow-band filters. The three upper panels represent results for the 
dynamic model Dl with mass loss, the three panels at the bottom are for the model Nl without mass loss. Triangles are for UD- 
radius computed taking into account only continuum opacity. For the continuum the same symbol is used for the different methods 
because there is no appreciable difference between the results. Crosses, stars and diamonds represent the different methods (one 
point, two points, least-square, respectively) used for fitting the visibility profiles of models with UDs. 



we can observe the periodic movement of the stellar interior, 
which is pronounced to a different degree for different models of 
TableQ] The periodic variation of the UD-radius in all the models 
without mass loss reflects the fact that the movement of the com- 
plete atmosphere follow mostly the pulsation. The time scales 
introduced by dust-formation in the atmosphere of the models 
with mass loss may deform the sinusoidal pattern, a behaviour 
that can be explained with outer layers which are not moving 
parallel to the inne r ones. This is well represented in Fig. 2 of 
iHofner et alj (f2003). but one should keep in mind that the move- 
ment of theoretical mass shells should not be directly compared 
to the changes in the atmospheric radius-density structure that 
we observe by interferometric measurements. 

From the lower three panels of Fig. [7] (no mass loss) it is 
obvious that there is not much difference between the UD-radius 
and UD-continuum radius for the model Nl, especially for the 
1.53 /urn and 9.975 //m filter close to minimum luminosity. We 
find the largest difference for the maximun luminosity (<p = 0) 
where the ratio between UD-radius and UD-continuum radius is 
1.09 for the filter 3.175 /an and 1.04 for the other two filters. 
For the model Dl (upper panels) the value of the same ratio is 
3.03 for the filter 1.53 fim, 4.29 for 3.175 fmx and 4.12 for 9.975 
/mi, much larger than the values found for the models without 
mass-loss. 

The three different methods of fitting give mostly the same 
result in the case of the model without mass loss, while in the 
case of models with mass loss, the UD-radius is smaller for the 



first 2 methods (one- and two-points fit) than for the least square 
result. This last result can be explained going back to our initial 
classification of profiles (Sect. [3]). The first and the third class 
of profiles belong to the "uncommon" profiles of S03, not ap- 
proximable with UD functions. In particular the third class of 
profiles (9.975 fim) is the one that mostly deviates from a UD 
(see Fig.[3}- 

Even increasing the number of fitting points, the resulting 
UD-radius cannot be considered reliable due to the particular 
shape of the profiles. From Fig. [7] it is also clear that for the 
narrow-band filters it is not possible to define a single, phase 
independent scaling factor that gives the possibility to convert 
from UD-radii to UD-continuum radii. 

We also computed the equivalent UD-radii and UD- 
continuum radii for the broad-band filters. The results are shown 
in Fig. [8] The broad-band has the effect of smoothing the vis- 
ibility profile and this can be noted in the UD-radii behaviour 
versus phase by the reduced amplitude of the radius-variation. 
As for the narrow-band also in the case of broad-band filters 
the UD-radius versus phase for the set of models D shows a be- 
haviour not synchronized with the pulsation. The ratio between 
UD-radius and UD-continuum radius is not constant with phase 
and/or with wavelength. For the phase zero it corresponds to 2.08 
in J, 2.83 in H, 2. 14 in K and 2.54 in U band. In the lower pan- 
els of Fig. [8] where the model Nl is represented, the resulting 
UD-radius is always closer to the UD-continuum radius. The ra- 
tio between UD-radius and UD-continuum radius for the phase 
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Fig. 8. UD-Radius versus phase for JHKL' broad-band filters. The four panels show the results for the model Dl with mass loss.Plot symbols are 
the same as in Fig. [7] 



zero is almost constant: 1.03 for J and H,\.0A for K and U . 
Again, the amplitudes of variation of the radius for the N mod- 
els are smaller than for the narrow-band filters, and the variation 
follows the near-sinusoidal behaviour of the inner layers. 

For the broad-band filters the difference in the fitting meth- 
ods confirms the trend found for the narrow-band case. For the 
models with mass loss the one-point and two-points fit meth- 
ods give smaller UD-radii compared with the one resulting from 
the least square method. The UD-radii, resulting from different 
methods, agree closely in the case of models without mass loss 
and the same result is obtained also for all the UD-continuum 
measurements. 

The dependence of UD-radius on the pulsation cycle is not 
covered by this contribution, it will be topic of a future contribu- 
tion. 

5.4. Comparison with M-type stars 

Compared to the M-type stars we emphasise the lack of a spec- 
tral window for measuring a pure continuum radius. A model- 
dependent definition of t he continuum radius i s thus needed for 
the cool Carbon stars. In llreland et al.l (l2004al) the phase varia- 
tion of filter radii in units of the parent star radius /?J3 is rep- 
resented for different near-continuum narrow and broad-band 
filters. They demonstrate how measurements of radii in near- 
continuum band-passes are affected by molecular contamina- 



3 R p is defined as the Rosseland radius that "the Mira variable would 
assume if the pulsation stoppe d and the stratification of the star becomes 
static" (Jacob & Scholz 2002). 



tions that mask the geometrical pulsation. This contamination 
is negligible only close to the visual maximum. Other suitable 
band-passes for measuring the pulsation of the continuum lay- 
ers for M-type stars are the broad-band filters J and H, or sub- 
bandpasses of H and K. 



To carry out a consistent comparison between the UD-radii 
co mputed in this work f or C-stars, and the resulting UD-radii 
of Jacob & S chord d2002l) for M-type stars is difficult because of 
the intrinsc difference in the modelling approach, and because of 
the different stellar parameters. But we can say as general state- 
ment that the ratios of UD-radii/UD-continuum that we obtain 
are systematically larger for our models with mass loss than the 
one pre sented for M-type stars in Figs. 4 and 5 of llreland et al.l 
(2004a). More comparable are the resulting ratios for the mod- 
els without mass loss. The largest difference is due to the strong 
molecular (first of all C2H2) and dust opacities that characterize 
the C-star spectra. The UD-radius increases with wavelength and 
in the case of Carbon-rich model atmosphere s with mass loss 
this in crease is larger than the one found by iJacob & Scholzl 
(|2002|) for M-type stars. In particular, the ratio between UD- 
radius in J and UD-radius in L' for our model Dl at phase zero 
is 2.07, while for the M - type m odels presented in Fig. 4 and 
Fig. 5 o f IJacob & Scholzl J2002) is 1 .4 in the case of largest dis- 
crepancy (phase 1.5). The same ratio for our Nl model is 1.03, 
meaning that this model without mass loss shows an increase of 
the UD-radius with wavelength lower than the one of the M-type 
stars. 
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Fig. 9. UD-Radius versus phase for JHKL' broad-band filters. The four panels show the results for the model Nl without mass loss. Plot symbols 
the same as in Fig. [7] 



6. Conclusions 

In this paper we present, for the first time, a theoretical study on 
the intensity and visibility profiles computed for a set of dynamic 
model atmospheres of C-rich AGB stars. The main results of our 
investigation can be summarised as follows: 

- The profiles computed in the narrow-band filters for our set 
of models can be divided in 3 morphological classes mainly 
according to wavelength region, with the exception of the 
filter located at 3.175 /im in the middle of the strong ab- 
sorption feature characteristic for C-star spectra. This filter 
shows a behaviour more closely related to the one of the 
narrow -band filters in the region around 12 /mi than to the 
surrounding wavelength range. We have chosen three filters 
to represent the morphology of intensity and visibility pro- 
files: 1.53,9.975 and 3.175 /im respectively for the first, sec- 
ond and third class. The intensity profile of the first class 
is characterized by a limb darkened disc followed by a tail- 
shape with a bump and a peak. The 3.175-like profiles show 
an extended plateau due to strong molecular opacity. The 
second class is intermediate between the other two. 

- The morphology of the profiles (narrow and broad-band) of 
models with mass loss is very far from the UD shape. Models 
without mass loss, and all the theoretical continuum profiles 
can be quite well approximated by a UD. 

- Bandwidth smearing affects the resulting broad-band profiles 
with an error of about 70% for certain baselines. For models 
with mass loss the difference is notable not only close to the 
first zero but also on the flank of the first lobe of the visi- 
bility profile. For a comparison of models with observations 



in broad-band it is important to take this effect properly into 
account. 

- C-star UD-radii are wavelength-dependent. The dependency 
is stronger in the case of models with mass loss, in this case 
the UD-radii show also large phase dependence. Around 
3 //m and in the N-band the star appears more extended due 
to C2H2 opacity. The UD-continuum radius is mostly inde- 
pendent of the wavelength and phase. 

- The extensions of the models increase for higher mass loss 
rate and lower temperature. 

- Models with mass loss show a complicated behaviour of 
the UD-radius versus phase compared to models without 
mass loss. This behaviour does not follow the pulsation of 
the interior. 

- It is shown that using only one or two points of visibility 
to determine the UD-radii of objects characterised by strong 
pulsation and mass loss, the radius of the star is smaller than 
the one measured using more points (least square method). 
In the case of models without mass loss the different methods 
of fitting give the same UD-radii. 

- The UD-radius is very close to the UD-continuum radius in 
the case of models without mass loss. 

- The trend of the UD-radii versus phase in broad-band is 
smoother than in the case of the narrow-band filters. 

- Contrary to O-rich stars, for C-rich stars no spectral win- 
dows for observing the continuum radius are available in the 
case of C-stars and some assumptions based on models are 
necessary. A useful replacement for a measured continuum 
radius is the UD-continuum radius defined in Sect. [5] which 
is basically independent of wavelength and pulsation phase. 
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- The difference in the UD-radii between the models should be 
detectable observationally, in particular in the L-band where 
MATISSE (lLopez et alj |2006) will provide new opportunity. 

- The radius computed with the UD function has to be cond- 
sidered only as a first guess for the real size of the star. The 
intensity and the visibility profiles of a C-star, especially in 
the case of models producing mass loss, are very far from 
being UD-like. 
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